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The linear dispersion relation for surlace perturbations, as derived by Levine et al, Phys. Rev. 
B 75, 205312 (2007) is extended to include a smooth surface energy anisotropy function with a 
variable anisotropy strength (from weak to strong, such that sharp corners and slightly curved facets 
occur on the corresponding Wulff shape). Through detailed parametric studies it is shown that a 
combination of a wetting interaction and strong anisotropy, and even a wetting interaction alone 
results in complicated linear stability characteristics of strained and unstrained solid films. 
PACS: 68. 55. J, Morphology of films; 81.15.Aa, Theory and models of film growth; 81.16.Dn, Self- 
assembly. 
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I. INTRODUCTION 



Studies of the morphological instabilities of a thin solid 
films arc a first step towards understanding complex phe- 
nomena such as the formation of a three-dimensional 
nanoscale islands in strained alloy heteroepitaxy. Such 
studies became common after the pioneering works of 
Asaro and Tiller flj, Grinfeld and Srolovitz [|. The 
classical Asaro-Tiller-Grinfeld instability is one of an uni- 
axially stressed solid film on a rigid infinite substrate. Its 
variants for the single-component and alloy films on the 
rigid as well as on the deformable substrates have been 
studied and this research continues. Reviews of works on 
the single-component films have been published, see for 
instance Ref. [4|. 

Film-substrate wetting interaction is a relatively new 
concept in the field of research on morphological insta- 
bility and evolution. When surface slopes are not very 
large, this additional wetting energy can be considered a 
function of the local film thickness h, but not the slopes 
of h 5-7]. In Refs. [sl-fioi [12} and others it has been 
shown that wetting interaction damps long-wave insta- 
bility modes in a certain range of film thickness, thus 
changing the instability spectrum from long-wave type 
to short-wave type. The latter mode of instability is 
more relevant to the process of formation of island ar- 
rays [13]. In Ref. Q it is recognized that in the pres- 
ence of wetting interaction, the boundary conditions that 
describe stress balance at the film free surface and at 
the film-substrate interface must be augmented by wet- 
ting stress terms - that are proportional to the rate of 
change of the surface energy with h. (Wetting stress 
is called conjoining (or disjoining) pressure when study- 
ing dynamics of thin liquid films on substrates (l^ [Ts} . 
This pressure is partially responsible for so-called spin- 
odal instability, which typically leads to film dewetting 
(rupture); for discussions of spinodal instability, see for 
instance Ref. [1^1 and references therein.) Wetting stress 
and lattice-mismatch stress have different origin, and the 



former may be present even when the latter is absent. 
For wetting (non-wetting) films, the solution of the elas- 
tic free-boundary problem with boundary conditions that 
include wetting stress terms, results in additional desta- 
bilizing (stabilizing) contributions in the dispersion rela- 
tion. Some stability characteristics have been analyzed 
in Ref. Q within the framework of longwave approxima- 
tion, where in addition the surface energy is assumed 
isotropic. This communication extends that work by 
adding strong anisotropy and considering not only wet- 
ting films, but also non- wetting films. Here we recognize 
that the film thickness and the wetting length are two 
independent characteristic lengths, i.e. the former length 
is determined by film deposition, while the latter one is 
determined by the molecular structure and properties of 
the film-substrate interface. Since wetting length may 
be, and normally is, less than the deposited film thick- 
ness, the perturbation wavelengths may be comparable to 
the film thickness but still much larger than the wetting 
length. In this case the long wavelength approximation 
may hold with respect to the wetting length, but not 
with respects to the film thickness. In Sec. IV we show 
that this approach reveals linear stability features that, 
we believe, went unnoticed in prior publications. 



II. PROBLEM STATEMENT 

Following Refs. d, [i3|, we consider a dislocation- free, 
one-dimensional, single- crystal, epitaxially strained thin 
solid film in a wetting interaction with a solid, semi- 
infinite elastic substrate. The film surface z = h{x, t) 
evolves due to surface diffusion. This evolution is de- 
scribed by 
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where f2 is the atomic volume, D is the adatoms difFu- 
sivity, N is the adatoms surface density, kT is the Boltz- 
mann factor, and Ai is the surface chemical potential 
[Is'] . The latter has contributions from the elastic energy 
in the film, the anisotropic surface energy, and a wetting 
interaction: 



cost 



,(2) 



where 9 is the angle that the unit surface normal makes 
with the reference crystalline direction, say [01] (cho- 
sen along the z-axis, which is normal to the substrate), 
7(/i, 9) is the height- and orientation-dependent surface 
energy, k is the curvature of the surface, S is the ar- 
clength along the surface and is surface slope (note, 

d/dS = {cos9)d/dx = {l + hiy^^^ d/dx). The term 
proportional to the small positive parameter S is the reg- 
ularization that is required in view of ill-posedness of Eq. 
^ for strong anisotropy, that is when > l/(m^ — 1) in 



Eq. (g]) below |19l - |22l |. Note also that the mixed deriva- 
tive term in Eq. ([2]) is nonlinear and thus it has no impact 
on linear stability. 

In this communication we consider the two-layer expo- 
nential model for the surface energy 0| : 

7(/i, 9) = jti9) + (7s - Itie)) exp {~h/£), (3) 

where 7s is the surface energy of the substrate when there 
is no film, £ is characteristic wetting length, and 7t(0) is 
the anisotropic surface energy of a thick film: 



7t(6l) =7o(l + 



e-v cos mt 



> 0. 



(4) 



Here 70 is the mean surface energy, is the strength 
of anisotropy and m is the integer parameter specifying 
anisotropy type (i.e., four-fold, six-fold, etc.) By com- 
parison with experimental and ab initio computational 
studies the two-layer exponential model has been shown 
the most accurate to-date [2^ In the absence of 

anisotropy, 7t = 70 = const., and S vanishes. S is taken 
zero also at weak anisotropy, < l/(m^ ~ !)■ In real- 
ity, the maximum of 7t(0) might occur a.t 9 — /3, where 
the non-zero angle /3 is a misorientation from the refer- 
ence direction. Without significant loss of generality we 
assume /3 = in Eq. (g]). 

The expression for the elastic energy £{h) in Eq. ([2]) 
is derived in Ref. [l7| without accounting for wetting 
interaction. Wetting interaction is considered in several 
papers, including Refs. [i-[ll [H [H [13 (in Ref. [13 the 
wetting effect arises not from the dependence of surface 
energy on thickness, but from the thickness-dependent 
elastic energy, which cannot be calculated from linear 
elasticity theory). To this end, by combining expressions 
derived by us in Refs. [1, [11] we state the dimensionless 
linear growth rate in the longwave limit, kho 1 (where 



k is the perturbation wavenumber and ho is the uniform 
thickness of unperturbed planar film): 

uj{ho, k, /i, e, Cj) — Ae^ (/i -|- Aihok) k^ — 
Be {ho -l) + ho {Biho - Ai) k] fc^aexp {-ho) + 
F [{A - (G + A) cxp {-ho j) k'^ - Afc*^- 
ak^ {exp {—ho) — -620 exp {—2ho)}\ .(5) 

I has been chosen as the length scale and ^^/Z? as the 
time scale. Here e is the misfit strain in the film, and 
/i = piij I \Xs is the ratio of the film shear modulus to the 
substrate shear modulus. Other parameters are: 



A = 



8iVf72as(H-i^/)V/ 



kTla) 



B = 



4A^r227oCks (1 + Vf) Vf 



kTP 



(6) 



^1 = 



2Ci 



= B, = A2^L, (7) 
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Ci — a-f + afl^sfi — a^lJ,'^ , C2 — 4:a f + 3a f Pgfi — ia^ii^, 

(8) 

= 2 (1 , = l-2i/y(,), (9) 
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A = (771^-1)6^-1, G = 7s/7o, 



G-1- 



(10) 
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In Eqs. (|6])-((9]) v is Poisson's ratio. Note coupling of 
wetting interaction and misfit strain through the term 
proportional to e (the second line of Eq. ([5])). This term, 
responsible for breaking symmetry between compressive 
and tensile stress states, drops out of the growth rate 
in the absence of wetting interactions, ho 00. In the 
square brackets of this term, — /i and —Aihok are the 
contributions from the wetting stress; see also Refs. [§, 
[25} . Another contribution from the wetting stress is the 
term proportional to exp {—2ho) in the square brackets 
in the last line of the equation. 

Our goal is to elucidate the roles of anisotropy, wetting 
interaction and wetting stress and to characterize film 
stability in the space of dimensionless parameters ho, k, 
fi, e and e^. Other material parameters will be fixed to 
their most characteristic values. We choose the following 
values: D = 1.5 x 10"^ cmVs, N = 10^^ cm'^, ^ = 
2 X 10^23 ^j^3^ f^r^ ^ X 10^13 gj.g^ ^ 2 X 10^ 

erg/cm^, i/j = 0.198, = 0.217, Hf = 10^^ erg/cm^, 
,5 = 5 X 10^12 gj.g, and £ = 3 X 10"^ cm. The value of 
the characteristic wetting length is of the order of 1 ML 
thickness for InAs or Ge film [1^. We assume strong 



anisotropy, i.e. > 1/ (rn^ — 1) 



and thus A > 



0. For strained films considered in Sec. IV, we choose 
777 = 32 as the most characteristic value [l^. However, 
as far as the effect of anisotropy on linear stability is of 
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interest, similar results are obtained for other common 
values such as to = 4 or to = 6. That is, choosing larger 
TO simply means that smaller values of e-y are above the 

critical value e^f . Wetting films require a > [2^ , thus 
we choose 7s — 2jq, ei^"^ < e-y < 1. For analysis of non- 
wetting films (a < 0), we choose 7^ = 7o/2, e-y > e'^K It 
is clear that wetting stress terms (pointed out above) are 
destabilizing (stabilizing) in wetting (non-wetting) films. 



III. FILMS WITH WETTING INTERACTION 
AND ZERO MISFIT STRAIN AND WETTING 
STRESS 

When misfit strain and wetting stress are not present, 
Eq. ([5]) reduces to: 

uiho, k, e^) = F [(A - (G + A) exp (-/iq)) k^ 

-Ak^ -ak^exp{-ho)] . (1) 

First, we consider wetting films. It follows that the 
perturbations with the wavenumbers larger than fcj, — 
K/ A cannot destabilize a film of any thickness. (Here, 
kc is not the customary cut-off wavenumber, but is deter- 
mined from the condition a; < for any /ig.) However, in 
the opposite case k < kc only the films of thickness that 



is less than the critical, h 



(ci) 
' 



ho < h^o'^ = -In- 



are stable: 



Afc2 _ Ak^ 
, + {G + A)k^ 



(2) 



With A = 25/9 corresponding to the material parameters 
stated above, m — A, and = 0.1, we obtain kc — 0.42. 
Taking typical k ^ O.lfcc in Eq. ^ gives h'-g'^ = 6.94, 
which translates to the dimensional value of 7 ML. Fig. [T] 
shows the contour plot of h'^^\k, e^). It can be seen that 

stronger anisotropy decreases h^^K We notice also that 
strong anisotropy destabilizes (that is, the contribution 
proportional to k'^ in the square bracket is positive) only 
relatively thick films, such that 



ho > h, 



ic2) 



'In- 



A 



G + A 



(3) 



For the chosen values, h^^'^'' = 1.6 ML. Such threshold- 
type influence of strong anisotropy is distinctly different 
from the simplified model in which wetting interaction 
is absent. The latter model can be obtained by taking 
ho ^ 00 va Eq. ([1]), and thus this equation becomes 
a;(A;,e^) — F (^Ak'^ — Ak^) , from which it is clear that 
strong anisotropy has destabilizing influence on a fllm 
of arbitrary thickness. These findings to some extent 
echo Refs. [23, 27], where the existence of the critical 
perturbation amplitude that is necessary to destabilize a 
film in the presence of a cusp in the surface energy 7(0) 
(which is the case below the roughening temperature), 
has been demonstrated. Thus if a film is thin, critical 
amplitude may be unattainable and the film will not be 




FIG. 1. Contour plot of the critical thickness /iq'^^' for strong 
anisotropy, > ei,"^' — 1/15. 



destabilized. However note that models of Refs. [23|, |27[ 
do not allow staightforward separation of the effects of 
surface energy and mismatch stress, and thus our results 
can't be easily compared to these papers. 

Next, we consider non-wetting films. One example 
of such material system may be the energetically-driven 
dewetting of silicon-on- insulator [2^, [2^ . Repeating the 
analysis and referring to the critical values shown above, 
it follows that film of any thickness is stable with re- 
spect to perturbations with wavenumbers larger than 



max 



(kc, fc^"^) , where fci"' = ^-a/(G-F A). If fc^ < fc < 

kc"\ then film is stable if Hq > h^^^^ and unstable other- 
wise. If fee"'' < k < kc, then the film is stable if hg < h^^^^ 
and unstable otherwise. Finally, li k < min (kc,ki"^^, 
then the film of any thickness is unstable. With G 



and e-y = 0.1, k, 



(u) 



0.5 

0.77 > kc, and therefore the second 
possibility, fci"^ < k < kc, must be dismissed. Typically, 
the first scenario (kc < k < ki^^) holds, and thus there is 
a critical thickness below which the film is unstable [TO] . 

Results similar to shown above for wetting and 
non-wetting films can be obtained (numerically) with 
non-zero wetting stress, since the negative exponent 
exp {—2ho) decays fast compared to the terms in Eq. ([T]) 
that are proportional to exp (—ho). 



IV. WETTING FILMS WITH NON-ZERO 
MISFIT STRAIN AND WETTING STRESS 

The situation presented in this section is common for 
Stranski-Krastanov growth of epitaxial thin films. 

As we pointed out in Sec. II, in the presence of mis- 
fit strain and wetting interaction, Eq. ([5]) contains the 
term that is proportional to the first power of misfit 
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strain. Whether this term is destabihzing or stabihzing 
(for, say, e > 0) depends on the sign of the expression 
/(fc, /i, /lo) = n{ho — l) + ho (Biho — Ai) k. Only for suf- 
ficiently small k and large this is positive. Then the 
second term in Eq. ([5]) is stabilizing as shown in Fig. [2] 
Increasing ho makes the domain of stabilization smaller. 
As fi is in the range 0.5 - 1.0 for a typical heteroepitax- 
ial semiconductor system, the coupling of misfit strain 
and wetting interaction has stabilizing effect on an ul- 
trathin film of thickness of the order of several wetting 
lengths, for longwave perturbations. Note that the stan- 
dard e^-term is always destabilizing for all perturbation 
wavelengths [l7| . 




0.20 



FIG. 2. Zero level curve of f(k,jj,,ho) — (fto — 1) -f 
ho {Biho - Ai) k. (a) ho = 2, (b) ho = 3. To the left of each 
curve, the symmetry-breaking term in the longwave growth 
rate ([5]) (second line) is stabilizing, to the right - destabilizing 
(when a, e > 0). 

In order to demonstrate some effects of arbitrary re- 
lation between wetting length, film thickness and the 
perturbation wavelength, in conjunction with strong 
anisotropy, we use next the full dimensional growth rate 
expression involving hyperbolic functions of kh, where 
k, h are now the dimensional wavenumber and mean 
thickness, respectively. (Available on request from au- 
thors.) (Eq. ^ emerges upon expansion of this growth 
rate in powers of small dimensionless parameter kh, re- 
taining the dominant terms of the expansion (longwave 
approximation), and non-dimensionalization.) As is Eq. 
[5l the full growth rate is quadratic in e, allowing one to 
explicitly determine the boundaries of neutral stability, 
w = 0, in the r — e or u — e planes. Here r and u (dimen- 
sionless) are defined hy h = r£, k — u{2n/£). 

Fig. [3]shows neutral stability curves, in u — e plane, for 
= and 0.01. (For value m = 32 used in this Section, 

e^'^'' = 0.001.) For all three values of a film thickness in 
the former (isotropic) case, and for the smallest value in 
the latter (strongly anisotropic) case, the film is destabi- 



lized by short-wavelength perturbations, u > 0.02, above 
some critical value of the misfit parameter e. Increasing 
film thickness in the isotropic case to values as large as 
50£ only makes the domain of stability to shrink. How- 
ever, for larger film thickness in the strongly anisotropic 
case (Fig. [3^0, d)) two stability domains emerge sepa- 
rated by the domain of instability. The splitting of a 
single domain into two domains occurs at r = 0.5. The 
size of stability domains decreases with increasing film 
thickness. Overall, the film is less stable with increasing 
anisotropy (as expected). Note that instability in Fig. 
[3fc,d) is present for some u even when misfit is zero. Re- 
sponsible for this is the combined destabilizing effect of 
anisotropy and wetting stress, which together overweigh 
the stabilizing effect of the wetting layer; also see Sec. 
III. Similar behaviour is observed for increasing /i while 
keeping thickness fixed. We also notice that only r = 0.1 
case can be (probably) captured by longwave approxima- 
tion, as kh = 2'Kru — 0.38 ~ 1 for r = 0.1, u = 0.6, and 
is even larger for other values of r in Fig. [31 In order to 
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FIG. 3. Neutral stability curves, (a); e^. = 0. r = 0.1 (solid), 
r = 1 (dash), r = 3 (dot), (b-d): = 0.01. (b) r = 0.1, 
(c) r = 0.518, (d) r = 1 (dashed), r = 3 (solid). Domains of 
surface stability (instability) are marked by S (U). 

characterize the horizontal spacing between two stability 
domains in Fig. [31[c,d), in Fig. 2] we plot the neutral 
stability curve corresponding to the level e = 0. It can 
be seen that for all reasonable r this spacing does not 
exceed 0.3. For comparison, the case of slightly larger 
anisotropy is also shown. 

To summarize, we considered all combinations of 
wetting interaction (through the exponential two-layer 
model), lattice- mismatch and wetting strains, and strong 
anisotropy. Our results demonstrate complicated linear 
stability of ultrathin films {h ^ 1 : b wetting lengths). 
In particular, we show that extremely thin {h ^ 1 : 2 
wetting lengths), unstressed wetting films are not desta- 
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0.00 0.05 0.10 0.15 0.20 0.25 0.30 
u 



FIG. 4. Neutral stability curves (see text), (a) e-, — 0.01, (b) 
e-, — 0.014. Film is stable below each curve. 



bilized by arbitrarily strong anisotropy. Anisotropic, 
stressed wetting films are destabilized by any level of 
mismatch stress, but only in the narrow range of per- 
turbation wavenumbers. Such films can remain sta- 
ble with respect to short- wavelength perturbations when 



they are very thin and at any reasonable mismatch 
stress level. Our final remark concerns two-dimensional 
surfaces and corresponding surface energy anisotropies 
"ftihxjhy) of the generic form Q (where contribution 
in the y-direction is additive, as is commonly assumed). 
We conjecture that, if the surface orientation (of a thick 
film) is still one of the high-symmetry crystallographic 
orientations, such as [001] or [111], then the efi'ect of 
such in-plane anisotropy is nonlinear and thus the lat- 
ter anisotropy will not affect the results. This can be 
qualitatively understood, for instance, by following the 
analysis leading to Eq. (15) in Ref. [Ij while accounting 
for the nonlinear nature of the mixed derivative term in 
Eq. 1^ and the form of Eq. Due to complexity 

of formulation and derivation, the exact proof is beyond 
the scope of this note. The results are also unchanged if 
the surface is two-dimensional but in-plane anisotropy is 
zero. 
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